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Abstract 

Approximations of the Sedov self-similar solution for a strong point explosion in a medium with the 
power-law density distribution p° cx r~ m are reviewed and their accuracy are analyzed. Taylor approxi- 
mation is extended to cases m 7^ 0. Two approximations of the solution are presented in the Lagrangian 
coordinates for spherical, cylindrical and plane geometry. These approximations may be used for the inves- 
tigation of the ionization structure of the adiabatic flow, i.e., inside adiabatic supernova remnants. 



1 Introduction 

Self-similar (Sedov ]l7| , |l8|] ) solutions for the strong point explosion in the uniform medium p° = const or in a 
medium with power-law density distribution 

~p°{f) = p°(0)r- m , (1) 

where f is the distance from the center of explosion, are widely used for modelling the adiabatic supernova 
remnants, solar flares and processes in active galactic nuclei. 

Sedov ([[l7| has obtained the exact solution solving the system of hydrodynamic differential equation 
on the base of dimensional methods. Independently, Taylor ([^9|) has solved the same task in the case of the 
uniform medium numerically and in analitical form approximately. The main Taylor's idea was to approximate 
the fluid velocity variation behind the shock front. 

Kahn (]Tl[|) have proposed the approximation of the Sedov solution in the uniform medium. His technic 
consists in approximation of mass distribution inside the shocked region. Using Kahn methodology, Cox & 
Franco (|(|) have built the approximation of the exact solution in the power-law medium ([!]) for m < 2. With 
the same technic, Cox & Anderson (Q) have presented the approximation for description of the shocked region 
and blastwave motion in uniform medium of finite pressure. 

Ostriker & McKee (p6[) basing on the virial theorem have given a number of approximations for the fluid 
characteristic variation as one- or two-power polinoms. 

Hnatyk (Q) have proposed to approximate firstly the relation between the Eulerian and Lagrangian coor- 
dinates of the flow elements. 

In present work, Taylor approximation is written for the medium with power-law density variation (Q). 
Using Hnatyk's approach, we develope also two approximations of the Sedov solution for power-law medium 
with to < 2 in Lagrangian coordinates that is useful for investigations of the nonequilibrium ionization processes 
in a shocked plasma, e.g., inside the adiabatic supernova remnants. One of the approximations presented here 
bases on the approximate hydrodynamic method for description of the nonspherical strong point explosion in 
the medium with arbitrary large-scale nonuniformity developed by Hnatyk & Petruk (JTo|]). Therefore, it may 
also be considered as additional test on this method. 



2 Sedov solution and its approximations 
2.1 Sedov solution 

If strong (P s / P° — > 00) point (R /R — * 0) explosion with finite energy E Q becomes in the point with coordinate 
f = in time t — 0, the blastwave creates and propagates with velocity D in the ambient medium with density 



p°(r) (P s and P° are pressure of the shocked gas and the gas of the ambient medium at the shock front position, 
R is the size of the body exploded, R is the radius of the blastwave). It is also assumed that injected mass is 
small and no energy lost from the shocked region during the motion. 

Such a task is described by a system of hydrodynamic differential equations. Sedov gives the analytical 
self-similar solution for description of the motion of shock front and the distribution of fluid parameters inside 
the shocked region for a strong point explosion in the uniform ambient medium and in the center of symmetry 
of a radially stratified medium ([!]). 

This solution shows that the strong blastwave in a medium with the power-law density distribution (|l|) 
moves with deceleration then m < N + 1 and accelerates then m > N + 1. If m > N + 1, both the mass inside 
any sphere, which containes the center of the symmetry, and kinetic energy equal infinity. We will consider 
m < N + 1 cases only. 

Radius R and velocity D of the strong blastwave in the medium (|l|) with m < N + 1 are (Sedov Jl8[): 

/ p \ l/(JV+3-m) 
R=( % J t 2/(N+3- m} (2) 

1 /2 

D(R) = — 4 (— ^ jr<"+ 1 ->/> ) (3) 



iV + 3-m Vcm p°(0) 

where N = 0, 1, 2 for plane, cylindrical and spherical wave, respectively, is a self-similar constant. 

Distributions of the fluid characteristics behind the shock front are self-similar, i.e., for any time t the 
density p, pressure P, fluid velocity u variations and coordinate a are 



p(f,t) = p s (t) 


■ P(r), 


(4) 


P(f,t) = p s (t) 


• P(r), 


(5) 


u(r,t) = u s (t) 


• u(r), 


(6) 


a(f, t) = R(t) ■ 


a(r), 


(7) 



where r = f/R(t), a is the original position of the fluid mass element and superscript "s" corresponds to values 
of the parameters at the shock front (Fig. [I]). 

Gas occupies whole shocked region (0 < f < R) when m < mi, 

mi= l + 3iV + (l-iV) 7 
7+1 

When m — *■ toi central pressure -P(O) — > 0. Shock waves in media with steep density gradients (m > mi) 
develop a cavity around the center of explosion. Such a cavity creates in the uniform medium (m = 0) when 
7 > 7i — (1+3N)/(N — 1). Sedov has also presented a solution for hollow blastwaves. Review of approximations 
for these cases is given by Ostriker & McKee (Jl6)). We do not consider m > mi in this paper. 
For m — mi (or 7 = 7* in the uniform medium) solution has very simple form: 

p(r) = r"" 1 , P(r) = r N+1 , 

(9) 



u(r) = r , a(r) 



= r (7+i)/(7-i) 



Singularities in the solution also appear with 7712 = (N + 1)(2 — 7) and 7713 = (27 + N — l)/ 7 then some 
exponents in the solution equal infinity. Similarity solutions for these cases are deduced by Korobejnikov & 
Rjazanov (@). For N = 2 and 7 = 5/3 mi = 2, m 2 = 1, m 3 = 13/5. 

Self-similar constant cca = ctA{N,j,m) in equations (§)-(|^) for R and D may be found from the energy 
balance equation with variations of density p, pressure P and mass velocity u inside the shocked region 

E _fp(r,tMr,t)^ Ndr fP(^ rNd ^ 



(7 J 2 J 7-I 







where a = 4tt for N = 2, cr = 2tt for TV = 1 and a = 2 for TV = or, generally, a = 2nN + (N - 1)(N - 2). If 
we proceed to normalized parameters using @-(^) and general shock front conditions 

Ps = ^ ft, Ps = -^TPsD 2 , Us = ~^—D (11) 

7 — 1 7 + 1 7+1 



we will obtain that E — (3a ■ MD 2 /2 with 

M = ap°(0)R N+1 - m /(N + 1 - m), (12) 



4(N+l-m) 

Pa = a — i (Ik + It) (13) 

7^ — 1 



constant shape-factor 
and constant integrals 

Ik = J P(r)u(r)*r n ' dr , J T = / P(/-j/-\/r . (14) 
o o 
Also we will have a self-similar constant 

2(7 

(N + 1 - m)(iV + 3 - m) 2 



l l 



°^ = , Ar | i , • &4 ■ (15) 



Simple formula gives <xa(N, 7, mi): 

2(7(7 + 1) 

c<A = 9 ■ (lb) 

(JV + 1)( 7 - 1)((N + 1)7 - N+ 1) 

The distributions (0)-(@) in the exact solution are parametric functions of an internal parameter. The 
expressions for the functions are complicated. These factors stimulate developing the approximations of the 
self-similar solution. 

2.2 Taylor approximation 

Basing on own numerical results, Taylor ( p9[ ) propose to approximate the velocity variation u(r) behind 
spherical (N = 2) shock front moving into the uniform medium (m = 0) as 

+ < (17) 

D 7 

where a and n are found to give exact values of u s , P s , p s and their first derivatives in respect to r. Substituting 
this approximation into the continuity equation and into the equation of state for perfect gas, the approximated 
distributions of the density and pressure obtain. Taylor do not give the dependence a(r), but it may be taken 
from the adiabaticity condition P(a)p(a)~ 7 = P(r)p(r) -7 and (|70|) - ([tT]) : 

a7 m-(iV+i) = p( r ) p ( r )-7 j (18) 

with approximations for P(r) and p(r). 

So, Taylor approximation for the variations of density p, pressure P, fluid velocity u and coordinate a are: 

*)-t-r"" (2±l.^iy, (1 „ 
Ps{t) V 7 7 / 

m .™.h±i-?iy, ,20) 

Ps(t) V 7 7 / 

aW = -| =r ,/«,-« (i±i_rri)-, (22) 

where n = (7 7 - l)/( 7 2 - 1), p = 2( 7 + 5)/(7 - 7), ? = (2 7 2 + 7 7 - 3)/(7 - 7), s = (7 + l)/(7 - 7). Self- 
similar constant = o^(2, 7, 0) goes with ([TJ]) and approximated profiles of p, P and w. Fig. [I] and table [| 
demonstrate accuracy of Taylor approximation in comparison with the exact solution. 
This approximation is extended to cases m =/= in section |3|. 
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Figure 1: a-c. Sedov solution and the accuracy of Taylor and Kahn approximations of the solution in the 
uniform medium: a exact Sedov solution, b relative differences of Taylor approximation, c relative differences 
of Kahn approximation. Lines: 1 - p(r), 2 - -P(r), 3 - u(r), 4 - a(r). 7 = 5/3. 



2.3 Kahn approximation 

Kahn apply his methodology to the strong spherical blastwave (N = 2) in uniform medium (m = 0) with 

7 = 5/3. It is proposed to approximate first the mass distribution 

r 

Kr) = = 3 / p{r)r2dr - (23) 


Sedov solution shows that P r (r) = near the centre (subscript "r" denotes a partial derivative in respect 
to r). This fact allows to find that fi r / // = 15/ (2r) at r — 0. On the base of the equation of motion, /i,-//i = 12, 
fi rr = 168 and (fx r /fi) r = 24 at r = I . Therefore ratio /i r //i is proposed to be approximated as 

^ = H + 9 - r \ (24) 

This formula satisfies all written boundary conditions at both ends. 
The mass distribution finds as integral from (pi]): 

M (r)=r 15 / 2 exp(A( r 8_i)). (25 ) 
Density distribution follows from (^) and (]25|): 

p(r) = /ir/3r 2 . (26) 



Adiabaticity condition gives pressure variation 



Velocity deduses from the mass conservation equation 

<r) = \r-^. (28) 

3 H r 

If present location of mass element a is r, then a(r) may be found from the condition of mass conservation 
/i(a) = jtt(r) and relation (fT3| ) ;u(a) = a 3 : 

a(r) =^W 1/3 ■ (29) 



The expressions for Kahn approximation are the same as (|3C|)-(|34|) with m = 0. The accuracy of this 
approximation are shown on the Fig. [l| 

2.4 Approximation of Cox &z Franco 

Appling Kahn's approximation technique, Cox & Franco (Q) obtain the approximation of the self-similar 
solution for an ambient medium with the power-law density distribution ([!]) with m < 2 for 7 = 5/3 and 
N = 2. Approximation of Cox & Franco are: 



p (r) = f | + 8 - 4m ) • r ( 9 - 5m )/ 2 



(30) 



5/3 



(31) 



x exp ( — (r 8 - 4m - 1' 

1 M(2-m) V ; 



1 i r 8— 4m 

a(r) = exp (^^y (r 8 ~ 4m - l)) , (33) 
Kr) = r 5(3- m)/2 exp ( r »-*» - 1)) . (34) 

Author's approximation for /3a is 

Pa = 1.125 ■ (0.22 + 0.52 • (3 - m)/3). (35) 
The accuracy of Cox & Franco approximation is shown on Fig. || and table [|. 

2.5 Approximations of Ostriker & McKee 

Ostriker & McKee (|L6)) in the frame of the virial theorem approach applied to spherical blastwave (N=2) 
in the power-law ambient medium (|l|) and time-dependent energy injection E a (t) oc t s , present a number of 
approximations for the self-similar solution. We consider further s = 0. 

Authors introduce the dimensionless moments of coordinate r and velocity u: 

l 

Kij =l„ I ' r*u(ry P {r)r 2 dr , (36) 



where = (7 + 1)(3 — m)/(7— 1), and consider three types of approximations for it(r) and p(r): linear velocity 
approximation (LVA) 

u(r) = r, p(r) = r ( 6 -^ +1)m >/(T- 1 J, (37) 

one-power aproximation (OPA) 

u(r)=r lu , p(r)=r l ", (38) 

and two-power aproximation (TPA) 



u(r) = a^A- 1 + (1 - a n )r z 



(39) 




Figure 2: a-c. Accuracy of Cox & Franco approximation of the self-similar solution in the power-law medium 
([l]): a relative differences of the approximation for m — —4, b relative differences for m = —2, c relative 
differences for m = 1. Lines are the same as on Fig. |l|. Fracture in the curves for p(r) and a(r) is due to very 
strong dependence of the relevant Sedov distributions on the internal parameter, which changes in these wide 
intervals of r on 10 -10 only. 



p(r) = apr 1 "- 1 + (1 - a p )r'p< 2 . (40) 

In such an approach the self-similar constant ola as well as exponents and l p may be expressed in terms 
of moments Kq 2 and K\\. Namely, under self-similarity aA — 2i:rj 2 Pa/ (3 — to), where r\ = 2/(5 — to) and factor 
/3a equals 

„ 2 2K 02 (3 7 - 5) + (5 -m)( 7 + l)^ii .... 
^3 (7 2 -D(7 + l) ■ (41) 

Exponents in OPA are 

2K m - Ku(l + K w ) _ 5K 2Q - 3 

(l-JfaoJ^n ' " 1-K.. 
Derivatives at shock front are used to obtain the moments. So 



1 

1 + s lj /l lJ 



K v = ■,,■/, . (43) 



where s.y = i + j in LVA and Sjj = i + j + — m)/2 in OPA. Using (43) may be written in a simple 
form in LVA: 

aA = 16 ^ ll7-5-m( 7 + l) (44) 

3(5 -m) 2 ( 7 a_i)(B 7 + i_ m ( 7 + i)) ' v ' 

Moments have more complicated form in TPA. In this approach the expression for u(r) coinsides with the 
approximation ( |2~l| ) of Taylor with n = 1 + 7 (mi — to)/( 7 — 1) that equals to Taylor's n at to = 0. So, TPA is 
extension of Taylor approximation of u(r) to to 7^ 0. Contrary to Taylor's approach to find p(r) and P(r) from 
the hydrodynamic equations, Ostriker & McKee find the density variation independently as TPA ( |Io| ) with 

7(7711 — m) 



y 10 - 7 - (7 + 2)to ' 

(45) 

3 — 7m 6 + (7 + l)(mi — 2to) 

ip,l = 7- j 'p,2 = : , 

7— 1 7 — 1 

where 7 > 1. For to = variation p(r) in TPA coinsides with the result of Gaffet (Qj) for case of uniform 
medium (Ostriker & McKee |16|]). Two-power velocity approximation is used to extend Taylor approximation 
to cases to 7^ in section ||. 
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Figure 3: Relative differences of Ostriker & McKee two-power approximation of the self-similar solution for the 
uniform medium. 7 = 5/3. Lines are the same as on Fig. [l]. 



Pressure distribution are also restored independently. It may be found in OPA as a linear pressure approx- 
imation (LPA) and for TPA in the frame of pressure-gradient approximation (PGA). 

Most of mass is concentrated near the shock front and distribution u(r) is close to a linear function of r. 
Therefore, as noted by Gaffet (@j), the right side of Euler equation 

9P(f,t) = 1 1 du(f,t) 

dM(r,t) 4vrr 2 dt ( ' 

is nearly a constant. LPA (Gaffet |Q, Ostriker & McKee |l6| ) use this feature assuming the pressure to be a 
linear function of the mass fraction fi(r) 

P(r) = P(0) + (P s */y M (r). (47) 

Logariphmic derivative of pressure at the shock front is P* = (dlnP/d\nr) s = (2'y 2 +7'j— 3— 7to(7+1))/(7 2 — 1). 
Mass in OPA is fi{r) = 3/"V<\ P(0) in LPA is (Gaffet §) 

P(0) = 1 + "° ; (48) 
oj(6 — m) 

where uf = u\R/D 2 = w((4 - 3w)(m - 3) + 2(1 -w)(4-2u- m))j2 (Hnatyk @), w = 2/(7 + 1). 

Such an approach (substitution with f~ 2 u\ instead of f~ 2 u t in (Eq)) was also used by Laumbach & Probstein 
((l5|) to develop the sector approximation. 

In PGA a power-law form for the pressure gradient 



dP{r) 



= P*^ 2 - 1 (49) 



dr 

is used to give two-power approximation for the pressure 

P{r) = P(0) + a p A 2 , (50) 



where a p = P*/Z p> 2 and 



P(0) = (7 + l)Vi-m) 



37 2 + 2O7 + 1 - (7 + 1)(37 + l)m ' 

(51) 

_ 3 7 2 + 2O7 + 1 - (7 + 1)(3 7 + l)m 
p ' 2 ~ 2( 7 2 - 1) 

Accuracy in determination of a a and P(0) in approximations of Ostriker & McKee is shown in table || and 
in revealing the flow parameters on Fig. |^. 

2.6 Cavaliere Sc Messina approximation of a a 

Cavaliere & Messina (||) with a simple technique approximate the equations for the radius and velocity of 
shock in the power-law medium (Q) and E a (t) oc t s . For s = his approximation gives 

'7-1 1 N + 1 - m 
7 2 - 1 V7+ 1 + 2 iV+ 1 



^ = -^T(^ + , 1 ^7^)- (52) 
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Figure 4: Accuracy of Hnatyk approximation of Sedov solution for the uniform medium. Lines: 1 - p(a), 2 - 
P(a), 3 - tt(o), 4 - r(a). 7 = 5/3. 

2.7 Approximate methods for an explosion in medium with arbitrary large-scale 
nonuniformity 

In this subsection we pointed out a number of approximate methods for description of a point explosion in 
arbitrary nonuniform medium. These methods may also be applicable for a medium with power-law density 
variation. Bisnovatyi-Kogan & Silich (Q) and Hnatyk (||) have given the reviews of these methods, their 
applications and accuracy. 

2.7.1 Thin-layer approximation 

Thin-layer approximation is firstly introduced by Chernyi (|Q) and used by Kompaneets (]l3|) and other 
authors to find analytical solutions for evolution of the shock front in a number of type of nonuniform media. 
It is assumed in this approach that all swept-up mass is concentrated in the infinitely thin layer just after shock 
front and the motion is stimulated with the hot gas inside the shocked region with uniform pressure distribution 
P(r) — 0.5 (excepting P s = 1). Layer of the gas moves with velocity u s . This method was developed to calculate 
anly the shock front dynamics and therefore does not allow to reveal the distribution of the fluid parameters 
behind the shock front. 

Thin-layer approximation gives for spherical blastwave in the uniform medium (Andriankin et al. 

aA = 75( 7 -l)( 7+ l)> ' (53) 

2.7.2 Sector approximation 

In the sector approximation, the characteristics of an one-dimentional flow find as decompositions into series 
about the shock front. 

Laumbach & Probstein ([jig]) have proposed the sector approximation applying it to spherical blastwaves 
in a plane-stratified exponential medium. Authors use Lagrangian coordinate a and propose to approximate 
pressure variation in the form equivalent to P(a) = 1 + -Pa(° — 1) (Hnatyk Q). Density variation is given by 
the adiabaticity condition and relation r = r(a) by continuity equation. Fluid velocity field is not determined. 
For shock radius and its velocity Laumbach & Probstein approximation yeilds in the uniform medium limit 

32tt(47 2 -7 + 3) 

aA = 225( 7 - 1)( 7+ 1)3 ■ (54) 

Gaffet ( @, §1 ) uses Lagrangian mass coordinates p, and finds pressure variation as a linear pressure approxi- 
mation P(/i) = l+P*(/i-l). Gaffet (@H) also propose to improve accuracy of the approximation, taking into 
account the second order coefficients in the series. Author calculates such coefficients in terms of Lagrangian 
mass coordinate p. Hnatyk (Q), considering different modifications of the sector approximation, presents the 
coefficients up to the second order in terms of a. 

2.7.3 Hnatyk approximation 

Hnatyk (|^]) introduces also the idea to aproximate firstly the relation r = r(a, t) between the Lagrangian a and 
Eulerian r coordinates of the gas element in each sector of shocked region. Density p, pressure P and velocity 
u variation behind the shock front are exactly deduced from this relation. Really, the continuity equation 

p°(a)a N da = p(f)f N dr (55) 



gives us the density distribution 

p{a) = im = im WD ) ' (56) 

the equation of adiabaticity 

P(a,t)=Kp(a,ty (57) 

yields the distribution of pressure 

P(M) = ( naW- 1 ( D(a)\ 2 ( p(a,t) 



and relation f = r(a, t) gives velocity 

u(a,t) 7 + 1 1 dr(a,t) 

u(a) = im=—D(R)^^- m 

Author propose to approximate r(a) as 

r(a) = a a exp (0(a - 1)) (60) 

with 

«-«) 2 -<a and p = rl a +rl-(r s a ) 2 . (61) 
Such an expression ensures the edge condition r(0) = 0, r s = 1 and values of the derivatives 

r s a = 1 - uj, (62) 
rJ =w(l-w)[3S + JV(2-w)-m] (63) 

where B = RR/R 2 , R = dR/dt is the shock velocity, m = —d\np°(R)/d\iiR, subscript "a" denotes a partial 
derivative in respect to a. 

This approximation is accurate near the shock front, but around the explosion site (for a < 0.1 or r < 0.4) 
characteristics do not restore correctly (Fig. ||). This approximation does not take into consideration any 
derivatives of r(a) near the center and the distributions of p(a), P(a), u(a) do not bind there, causing such a 



situation. This approximation is extended to the central region in subsection 4.5 



3 Extension of Taylor approximation to m^O 

In this section Taylor's technic is applied to the case of a medium with the power-law density distribution (|l|) 
with m < mi. Ostriker & McKee ( fl6|) give the coefficients in the approximation ( |l7| ) for u{r) in such a case. 
Approximated velocity variation is (|2l|) with n = P* — 2 = (77 — 1 — 7717(7 + l))/(7 2 ~ !)■ Substitution with 
(|l7) into the equations of continuity and state gives 

p r 3 + 07(71 + 2)r"~ 1 — 7717 



p (7 — l)r — orfr n 

P r _ a~i 2 (n + 2)r n - 1 
P (7 — l)r — o;7r n 
After integration, pressure variation will be expressed with ( pp| ) where 

27 2 + 77 - 3 - 7717(7 + 1) 



7 — 7 — 711(7 + 1) 



Density is 



where 



Eq. (E2|) gives a(r) with 



7 7 

2(7 + 5-771(7 + 1)) 
7 — 7 — 711(7 + 1) 

7 + 1 



(64) 
(65) 

(66) 



p(r) = r (3—7)/(7-D (l±i_l__) ( (67) 



(68) 



( lT • (69) 
7 — 7 — 711(7 + 1J 

Exponents 71 = 6 — 5m/2, q = (16 — 5m)/ (3(2 — m)), p — (5 — 2m)/ (2 — m) and s = 1/(2 — m) for 7 = 5/3. 
This extended Taylor approximation is compared with the exact solution on Fig. 
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Figure 5: a-c. Accuracy of the extended Taylor approximation of the self-similar solution in the power-law 
medium (|l|): a relative differences of the approximation for m = —4, b relative differences for m = —2, c relative 
differences for m = 1. Lines are the same as on Fig. [l]. 7 = 5/3. 

4 Approximations of the Sedov solution in Lagrangian coordinates 

In this section we present two analytical approximations of the self-similar solution for a medium with the 
power-law density distribution expressed in Lagrangian geometric coordinates a. 

4.1 Flow characteristic distributions 

Exact expressions for normalized density p and pressure P variations behind the shock front moving into the 
power-law medium (p|) follow from (|5^) and (|58|) : 

(70) 



p(a) = ^-a N - m -(r(a) N -r a (a)) 1 , 
P(a)= 07Tl) 7 ■* N( - < - 1) - 1 -(r(a) N -r a (a)y 



(71) 



Distribution of the fluid velocity it (a) may be found from (|59|). Due to f = rR time derivative df/dt = 
Rr t + Rr a a t + rD (subscript "£" denotes a partial derivative in respect to t). We have also that at = —aD/R 
and, in the self-similar case, r t — 0. So, 



u(a) = 



7+1 



r(a) - r a (a)a 



The distribution p(a) follows from the definition (|23|) and (|55|): 

M (a) = < W+1 >- TO . 

4.2 Self-similar constant 

Self-similar constant o,a(N,^, to) in equations for R and D (§)-(||) obtains from ( [l5|) and (|l3|): 

8 (7 



(72) 
(73) 



OtA 



7 2 -l (3 + iV-m) 5 



(/k + Jt) , 



with 



- [ (r(a) - r a (a)aj a 



N—m 



da 



1 

(7Ti)7( r(a)JVra(a) )" V(7 " 1) "" a - 



(74) 



(75) 



Tabic 1: -P(O) calculated according to the self-similar solution and C for a strong point explosion in the 
power-law medium. 



N m P(0) C A 







7 = 7/5 


7 = 5/3 


7 = 7/5 


7 = 5/3 








0.3900 


0.3532 


1.1429 


1.1670 


1 





0.3729 


0.3215 


1.0863 


1.1112 


2 





0.3655 


0.3062 


1.0618 


1.0833 




-4 


0.4268 


0.3954 


1.0233 


1.0293 




-3 


0.4193 


0.3848 


1.0276 


1.0350 




-2 


0.4088 


0.3696 


1.0339 


1.0433 




-1 


0.3928 


0.3463 


1.0438 


1.0570 




1 


0.3087 


0.2217 


1.1054 


1.1556 




2 


0.1273 


0.0000 


1.3648 


1.0000 



4.3 Factor C and exponent x 

In Sedov self-similar solution, if r — > then the dependence r(a) is 

r = C-a x . 



(76) 



For m < mi, if we substitute (76) into (|7l| ) we obtain the connection between the factor C and normalized 
central pressure P(0): 

1/(JV+1) 



C 



7 



7 + 



r P(0) 



-1/7 



(77) 



We have to put x — (7 — l)/7 during this transformation in order to satisfy condition P(0) 7^ 0. In the case 
m = mi the exact solution (0) gives x and C = 1. General formula for exponent x is 



(7 — l)/7 for m < mi 
(7 — l)/(7 + 1) for m = mi 



(78) 



Analytical expressions for P(0) from self-similar solution are presented in Appendix 0. Calculated values 
of P(0) and C for a number of N, 7 and m are shown in table [j]. 



4.4 Derivatives at shock front 

Expressions for the derivatives r%, r s aa , r s aaa may be obtained with the technic of Gaffet (0) from the set of 
hydrodynamic equations for perfect gas and conditions on the shock front (see Hnatyk & Petruk (|l(i||) for 
details). Derivatives r a , r aa are given with (62)-(p3|) and 

7t aa =^(l-^)[3(7-5c)P 2 + 

+ [(-5w 2 + 4w + 8)N + (4w - ll)m] P+ 
+lo(2lu 2 -7lu + 6)N 2 + (lu 2 +lu - A)Nm- 
-w(2 - w)^V - (w - 2)m 2 + (2w - l)m+ 
+(2w- l)m' + (6w-4)q1, 



(79) 



where Q = R 2 R^/R 3 and m' = -dm/ din R. 

In the power-law medium (Q) rn' = 0. Taking into consideration the equations for the shock radius (||) and 
shock velocity (||) we may also write 



B = - 



N - m + 1 



(N - rn + 1)(N - m + 2) 



2 ' " 2 
Reduced expressions for the derivatives r s a , r s aa , r s aaa are shown in table I 



(80) 



Table 2: Derivatives of the relation between Lagrangian and Eulerian coordinates at shock front moving into 
the power-law medium (fil). 



Derivative 



7 = 7/5 



r- = i 
a 6 



(-27V + 3to-9) 



2 3 3 3 
5 



(-2N 2 + 9Nm + 1837V - 9m 2 - 270m + 675) 



7 = 5/3 



a 4 

r a aa = ^(-N+2m-6) 

rlaa = J> {Nm + 19N - 2m 2 - 36m + 94) 



4.5 Second order approximation 

So, to approximate the self-similar solution, we approximate the relation r — r(a) between Eulerian r and 
Lagrangian a coordinates of flow elements. Following to Hnatyk's approach (|6(]) and like to relation j3^), 
r = r(a) may be aproximated in the form 

r(a) = a x exp (a{a p - 1)) (81) 

with x given by ( |78| ) and 

(r£ - x) 2 _ r°„ + < - (r^) 2 



or, after substitution with (|62j)-(63), 

2(1 -lu-x) 2 



lo(1 - uj)(N + m - 1 - 27Vw) ' ^ 
(3 = a~ 1 (l-w - x). 

Such a second order approximation, besides r(0) = 0, r s — 1, r„, r^ a , gives (91nr/91na)° = x, and, contrary 
to Hnatyk approximation, extends description of a flow to the central region. 

Variations of p(a), P(a) and u(a) follow from (|70|)-(|72j). For case N — 2, 7 = 5/3 and m < 2 these relations 
give (3 = 5(2 — m)/8 and 



(85) 



p(„) = f | _ |o?) - (8-!")/5 



lO—i 

«" o - 

3 
13 
CO 



-10 - 



CD 



-20 





Figure 6: a-c. Accuracy of the second order approximation of the self-similar solution in the power-law medium 
(0) for 7 = 5/3 and N = 2: a relative differences for m — 0, b relative differences for m = —2, c relative 
differences for m = 1. Lines are the same as on Fig. ||. 



Approximation (p^)-(p7[) may be considered as an inversion of Cox & Franko approximation (^0[)-(p3[). Unfor- 
tunately, accuracy of presented formulae is lower (Fig. ||, table ||) . 

4.6 Third order approximation 

In order to improve accuracy, we postulate the approximation r — r(a) to give exact values of two additional 
derivatives: third order r s aaa and (dr/d(a x ))° — C. Consideration of r^ aa is equivalent to consideration of the 
second order derivatives p s aa , P^ a , u s aa in expansion of relevant characteristics into the series near the shock 
front. This approximation is the same as used in the approximate hydrodynamical method for modelling the 
asymmetrical strong point explosion in the medium with a large-scale density nonuniformity (Hnatyk & Petruk 
[p)[). Contrary to the method, we take here that both the self-similar constant a a and factor C are different 
for different to. 

Namely, if at time t the shock position is R(t), we approximate a connection r = r(a) as follows 



r(a) 



where £ = 1 — a. Coefficients a, (3, S and exponent x are choosen from the condition that the partial derivatives 
at the shock front (a = 1) as well as (<91nr /<91na)° = x and (dr/d(a x ))° — C in the place of 



a > aa ' aaa 



explosion (a = 0) equal to their exact values: 



a = — r„ + x, 



In terms of a relation 



P = ~-{r s aa -2x-r s a + x(x + l)), 

S = g ' ( — r aaa + %X ■ T aa ~ 

-3x(l + x) ■ r s a + x(x + 1) [x + 2)) , 
<5 = C-(l + a + /3 + <r). 
i) and its first derivative are 

r(a) = a x (B - B x a + B 2 a 2 - B 3 a 3 + B 4 a 4 ), 



(89) 



(90) 



with 



r a (a) = a x-1 (A) - Ma + A 2 a 2 - A 3 a 3 + A 4 a 4 ), 



(91) 



Tabic 3: Self-similar constant ap i _{N^ r y,m) calculated with third order approximation of r(a) (88). 



N m aA 







J — I/O 


7 — 0/6 








1.0763 


0.6018 


1 
1 


n 
U 


u.yo4i 


U.OD44 


2 





0.8519 


0.4944 




-4 


0.2295 


0.1270 




-3 


0.2960 


0.1650 




-2 


0.3966 


0.2232 




-f 


0.5598 


0.3192 




1 


1.4631 


0.8722 




2 


3.3537 


1.8235 



Table 4: Coefficients in approximation (^0|). 7 = 5/3. 



N 


??? 


S 




B 2 


B 3 


Bi 








1.1670 


0.1333 


-0.1127 


-0.1074 


-0.02833 


1 





1.1112 


-0.01510 


-0.24655 


-0.1530 


-0.03276 


2 





1.0833 


-0.05189 


-0.2130 


-0.08708 


-0.009294 




-4 


1.0293 


0.2837 


0.9302 


0.9766 


0.3008 




-3 


1.0350 


0.1744 


0.6152 


0.6971 


0.2213 




-2 


1.0433 


0.07170 


0.3041 


0.4162 


0.1406 




-1 


1.0570 


-0.01334 


0.01359 


0.1452 


0.06128 




1 


1.1556 


0.08965 


-0.1753 


-0.1471 


-0.03778 




2 


1 















B = l + a + (3 + (; + 8 = C, A Q = xB Q , 

Bx = a + 2(3 + 3c: + 4(5, A 1 = (l+x)B 1 , 

B 2 = + 3c: + 6(5, A 2 = {2 + x)B 2 , 

B 3 = c: + 4(5, A 3 = (3 + x)B 3 , 

B 4 = S, Ai = {4 + x)B 4 . 

Distribution of p(a), P(a) and u(a) obtain from (|70|)-(|72]). Self-similar constant a a is given with (|74|). To 
simplify the procedure, numerical values of cxa(N, 7, m) in this appoximation are presented in table |[ Table ^ 
gives ready-calculated values of the coefficients in the approximation (|9^) for a number of cases. 

The accuracy of flow characteristic distributions in this approximation is high for uniform medium (Fig. [?]). 
Approximation coinsides with the exact solution (^|) for case m — mi. For other m 7^ 0, differences increase 
with increasing \m\ but maximal errors reveal in the region with low densities (Fig. |^). We compare also 
numerical values of a a and P(0) in this approximation with those from exact Sedov solution in table |^. a a in 
the approximation is close to the exact values and gives accurate shock radius R and velocity D. 



5 Conclusions 

In this paper, we review approximations of the self-similar solution for a strong point explosion in the power 
law medium p° oc r~ m and compare their accuracy with the exact Sedov solution of the problem. Different 
approaches found on the different basic approximations. Namely, Taylor (|l9)) and Ostriker & McKee (|T^|) 



I I I I I I I I I I I 
0.0 0.2 0.4 0.6 0.8 1.0 

a 

Figure 7: a-c. Accuracy of the third order approximation of the Sedov solution in the uniform medium (m = 0) 
for 7 = 5/3: a relative differences of the approximation for N = 0, b relative differences for N = 1, c relative 
differences for N = 2. Lines are the same as on Fig. |[ 

Table 5: Comparision of the self-similar constant aA and pressure P(0) calculated with: S - Sedov (Cffl) 
solution (Kestenboim et al. T - Taylor (JljJ) approximation; CF - approximation of Cox & Franco (||); 

LVA, OPA and TPA of Ostriker & McKee (Q); CM - approximation of Cavaliere & Messina (f|); TL - 
thin-layer (^3|) approximation; LP - approximation of Laumbach & Probstein (|l5|]); SOA - second order 
and TOA - third order ( |S8|) approximations. Uniform medium, 7 = 5/3 and N = 2. 



S T CF LVA OPA/LPA TPA/PGA CM TL LP SOA TOA 

a A 0.4936 0.4957 0.4930 0.5386 0.5027 0.4957 0.5655 0.5655 0.4398 0.4981 0.4944 

P(0) 0.3062 0.2855 0.3140 0.3333 0.3333 0.5000 0.3333 0.2507 0.3062 



approximate firstly the fluid velocity variation behind the shock front. Taylor used approximated u(r) sub- 
stituting it into the hydrodynamic equations to obtain full description of the flow. Contrary to this, Ostriker 
& McKee approximate p(r) and P(r) independently. Kahn (Q) technic, used also by Cox & Franco (@), 
consists in approximation of the fluid mass variation fi(r) and further usage of the system of hydrodynamic 
equations. Gaffet (0), Laumbach & Probstein (]lq]), Ostriker & McKee (Q) base their approaches on the 
approximation of or P(r). Thin layer approximation may also be included into this group. Hnatyk 

(J9|) take approximation of the connection between Eulerian and Lagrangian coordinates as basic relation. So, 
practically all possible approaches are used to have approximation for the self-similar solution. 

In this paper we apply Taylor's methodology to discribe a strong point explosion in the power-law medium, 
extending his approximation written for uniform medium, and write also two approximations expressed in 
Lagrangian geometric coordinates, approaching r(a) with different accuracy. 

Errors of all approximations are caused only by errors in the basic approximation. When the first approxi- 
mation has higher accuracy we have more accurate approximation for parameters of the shock and flow. 



Appendix: central pressure P(0) 

In this appendix, we give exact expression for P(0) in self-similar solution when m < mi (Sedov ]l8fl ) and 
when m = m 2 (Korobejnikov & Rjazanov 0). These relations complite the full set of formulae to build the 




0.0 0.2 0.4 0.6 0.8 1.0 



a 

Figure 8: a-c. Accuracy of the third order approximation for power-law medium, 7 = 5/3 and N = 2: a relative 
differences for m = —4, b relative differences for m = — 2, c relative differences for m = 1. Lines are the same 
as on Fig. ||. 



third order approximation of the Sedov solution for any 7, m < min(iV + 1, mi) and type of symmetry (plane, 
cylindrical or spherical blastwave). 

P(0) = for m = m x . 

In the case of m < mi and m ^ 1x12 

m =(iY(2±±Y(^r , (92) 



£1 



£2 



If m = m,2 then 



2 y \ 7 J \m — mi 
2(N+l) 



N + 3-m ' 

2(iV+l) 7 (iV+l-m) 



iV + 3-m (JV+l)(2-7)-m 



_ (N + 1 - m)(N + 3 - m) 
£3 = (JV + l)(2- 7 )-m + m ~ 2 ' 



e 4 



7 + 1 2 



(iV+l)(7-l) + 2 AT + 3-m 

7-1 



7(2 - m) + N - 1 



P(0)=( 5 ) (1±1) exp(-2 £ , , ( 93 ) 

2(iV+l) 



(JV + l)(7-l) + 2 ' 
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